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ABSTRACT 



Aims. We quantify the mixing of tiie measured cosmic-shear E- and B-modes caused by the lack of shear-correlation measurements on small 
and large scales, arising from a lack of close projected galaxy pairs and the finite field size, respectively. 

Methods. We calculate the aperture-mass statistics (M^ j^> and the E-/B-mode shear-correlation functions ^e.b± where small- and large-scale 
cutoffs are taken into account. We assess the deviation of the obtained E-mode to the true E-mode and the introduction of a spurious B-mode. 
Results. The measured aperture-mass dispersion is underestimated by more than 10% on scales smaller than 12 times the lower cutoff. For a 
precise measurement of the E- and B-modes at the percent level using a combination of ^e,b+ and ^e,b-, a field as large as 7 (2.4) degrees is 
necessary for ground-based (space-based) observations. 
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1. Introduction 

The observation of the correlation between the shapes and 
orientation of high-redshift galaxies has become a well- 
established method of studying the dark matter distribution 
on very large scales. Using theoretical predictions for the 
power spectrum of the projected cosmic density field, cos- 
mological parameters can be obtained from the measure- 
ment of cosmic shear, most notably the matter density Qm 
and the amplitude of densi ty fluctuations cth. Very recent re- 
sults i r iclude ground-based (iHoekstr a et alJ2 005UMassev et al.l 
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2005; Semboloni et al. 2003; Jarvis et al. 200(>) and space- 
based ( Rhodes et al. 2004: Hevmans et al. 2005) observations. 
Presently, cosmic-shear surveys are still limited by system- 
atic errors arising from the imperfect shape measurement of 
faint galaxies and a deficient PSF correction. Extensive stud- 
ies and compari sons between different data analysis methods 
are being made JHevmans et al.ll2006al) to find sources of any 
contamination of the cosmological signal. A common means 
to check for systematics in the data is the decomposition of the 
measured shear field (or power s pectrum) into th e gradient- and 



curl-parts (E- and B-modes), see lCrittenden et al.it2002!) . here 
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after C02, and lSchneider et al.1 ( l2002l) . hereafter S02. Since the 
image distortions caused by gravitational lensing are (to first 
order) curl-free, the presence of a B-mode is a distinct imprint 
of residuals not completely removed from the measurement. A 
B-mode of cosmological origin can be caused by the iiitrinsic 
correlation of galaxy orientations (e.g. ICrittenden et alJl2001t 
|jing|2002h . However, this effect is thou ght to be very small for 
reasonably deep surveys (e.g. Hirata e t alJl2004l) . and it can be 
removed from the shear signal using photometric redshifts or 
by modeling the intrinsic alignment signal that is distinct f rom 
the shear sig r ial ( Hevman s & Heavens 2003; King & Schn eideii 
l2002i l2003t lKing|l2005l) . A further non-negligible source of 
confusion could be the recently discovered intrinsic shape- 
shear correlation (jHir ata & Seliak 2004). Th is contamination 
also creates a B-mode JHevmans et alJ2006a) . and it is not yet 
clear how to correct for this effect. 

As we present in this paper, up to now the E- and B- 
mode decomposition cannot be performed directly from the 
shear correlation as measured from the data, but it involves an 
integral over the shear-correlation function ^-t to either arbi- 
trary small or infinitely large angular separations. However, the 
scales on which ^± can be measured are limited. On arc-second 
scales, the blending of closely projected galaxy pairs prohibits 
a reliable determination of the shape of those galaxies. On large 
angular scales, the measurement is restricted by the finite field- 
of-view. These limits cause a mixing of the E- and B-modes 
with the currently-used estimators, preventing their clear-cut 
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separation. Although this is a well-known fact, to our knowl- 
edge this mixing has not or only inadequately been taken into 
account to date. 

In this paper, we quantify the mixing of the E- and B-mode 
due to the small- and large-scale limits of the shear-correlation 
measurements. In Sect. |2 the aperture-mass dispersion and in 
Sect. |3] the E-/B shear-correlation functions are considered as 
measures of the E- and B-mode, respectively. Sectionl^briefly 
discusses the dependence on cosmology. In Sect. |5] we give a 
summary and offer ways of clearly determining the E- and the 
B-modes. 

2. Aperture-mass dispersion 

The aperture-mass statistics (IKaiser et al.l Il994t ISchneiden 
Il996l) is a smoothed function of the convergence field, 



Map(0) = idHUemX^), 



(1) 



where Ug is a compensated filter function, i.e. J d§&U0(&) - 
0. Equation Q can be written as an integral over the 
tangential shear jt, weighted with the function Qe{'&) = 
2/i?2 j^ d&' §' Ue{9') - Ue{&). Moreover, M^ is defined as the 
weighted integral over the cross-component of shear y^, 

M,p,AO) = j d-^Qe(m)yrA^). (2) 

The dispersion of Q dSchneider et alJl998l) can be written 
in terms of the convergence power spectrum f ^, 

2^J 



(Mt^m 



de£P^{{)u\ee), 



(3) 



where Uirf) is the Fourier transform of the filter function 
ui-d/O) = Q-Uei'd). The dispersion of M^ is non-zero only if 
a curl-part or B-mode is present in the shear field. This is not 
the case for a purely gravitational shear signal (to first order). 
Therefore, the measurement of {M\) is a test for systematic 
errors and/or intrinsic galaxy alignment. 

Two sets of filter functions have been extensively used for 
cosmic-shear measurements and studies. One set are polyno- 
mial functions with finite support dSchneider et al. 1998,) . 
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The 2PCF are defined in terms of the tangential and cross- 
component of shear, yt and jx, respectively, where these two 
components are measured with respect to the direction con- 
necting two points. The 2PCF can be written as functions of 
the power spectrum as follows. 
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where H is the Heaviside step function. The other set is related 
to Gaussians (C02), 
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The aperture-mass dispersion can be obtained by integrat- 
ing over the two-point correlation function (2PCF) of shear, ^-t 
(C02, S02), 
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^±W = (rtr,)('?)±<rxrx)W= dtEPMioAm- (7) 

Jo 
The functions T± in (|5) are given by 



r±W = — dtthA{xt)U\t); 

27T Jo 



(8) 



closed-formula expressions in the case of the above defined fil- 
ter functions can be found in S02 and C02. Expression (|6j is 
very useful since direct determination of the aperture-mass dis- 
persion from data is very inefficient due to boundary effects and 
masked regions in the images. The 2PCF, on the other hand, can 
be obtained easily from any survey topology. Note that, in the 
case of the polynomial filter, the integral in (|6j only extends to 
26, since then T±(x) - for x > 2. 

Projected close galaxy pairs in the data have to be excluded 
from the shear measurements. The shape of close pairs can- 
not be estimated reliably because of the blending of the in- 
dividual galaxy images. That means that the 2PCF cannot be 
measured for smaller angular separations than some cutoff sep- 
aration fljnin- For space-based observations, this cutoff may be 
as small as one arc second. Ground-based observations usually 
are limited by flmi,, ~ 5 arcsec. 

This inevitable cutoff at small angular scales leads to a 
biased estimate of the aperture-mass statistics. We introduce 
the quantity (M^ J,Q, flmin)), which is defined as in (|6}, but 
with the lower integration limit replaced by O^i^. Even in the 
absence of a B-mode in the convergence field, a O^i^ > 
will create a spurious B-mode signal. For small angular scales 
§, the first term in square brackets in (|6j is approximately 
a positive constant, since both ^+ and T+ are integrals over 
a positive function multiplied by Jq. The second term tends 
to zero because of J4 in both ^_ and T . This will result in 
a negative B-mode, {M\(0,0inm)} < 0, and an underestima- 
tion of the E-mode, {Mlp{0, 0mm)} < {MI^(_0)). In fact, since 
the r_-term is small for small scales, the approximate relation 
{Ml^{0, 0„,i„)) - {Ml{0, 0™i„)) ^ {M%{0)) holds for 0m,J0 « 1. 

Note that the same effects would arise if {M^ ^) were ob- 
tained by putting apertures on the data field. In this case, the 
estimator of the aperture-mass dispersion involves summation 
over all galaxy pairs in an aperture. Due to the cutoff, this sum 
would lack close pairs and the E-mode would be underesti- 
mated. 

Figures[2and|2]show the influence of the cutoff ^mi,, due to 
the absence of close pairs on the E- and B-modes; see Sect.|4] 
for the dependence on the cosmological model and source red- 
shift distribution. In Fig.|3]we quantify the suppression of the 
true E-mode. Only for aperture radii larger than the threshold 
0Q is the deviation smaller than indicated by the curves. For the 
polynomial filter, one roughly finds 0q - 120min for 10% accu- 
racy. For example, for ground-based observations with 0„^i„ = 5 
arc seconds one gets deviations of more than 10% for 5 1'. If 
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Fig. 1. The E-mode {ML{6)) and the leakage from the E-mode 
(M3p(0, flniin)) into the B-mode (Mj(0,0„in)) due to the small- 
scale cutoff e^in. The ratio <M^p(0, 0min)>/<A^ap(^)> is plotted in 
the small panels. 



a 1%-precision is aspired scales smaller than 3'7 have to be 
discarded. 

The Gaussian filter is less affected by the cutoff than the 
polynomial one, since it is broader and samples the 2PCF on 
larger angular scales for a given aperture radius. However, 
this advantage here turns into a disadvantage on large angu- 
lar scales. There, the E- and B-modes cannot be determined 
reliably due to the field boundary and the infinite support of the 
Gaussian filter functions. 

The definition of r+ (jS) shows that the mixing of the E- and 
B-modes due to a small-scale cutoff is inevitable. Regardless 
of the choice of the filter function, Ue, T+{x) tends to a positive 
value for x — > 0. The more rapidly T+ falls off, the smaller the 
bias and thus the shallower U and the broader Ue is. 

3. E- and B-mode shear-correlation function 

The E- and B-mode correlation functions ^e,b+ and ^e,b- were 
defined in S02 and C02. They can be obtained from the 2PCF 
in the following way, 

ft±(0) = \ [Ue) + €-(0) + ^u,l(0)] ; 
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2. The B- to E-mode ratio for various small-scale cutoffs 
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Fig. 3. On scales 9 larger than 6q, the ratio 
{Mlp(e,0mm))/{Mlp(e)) is lower than 1 - p, with hnes 
corresponding to p - 10%, 1%, and 0.1% plotted as a func- 
tion of the cutoff scale 0min- Thick curves correspond to the 
polynomial filter @, thin lines to the Gaussian functions (|5|i. 



ft±(0) = ^[^+(0)-^(0)+^u,l(0)]. 



(9) 



where the integral functions ^u and ^l correspond to the '-i-'- 
and '-'-cases, respectively, and are given by 
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Thus, in order to separate the E- from the B-modes one 
needs to know either ^_ up to very large or ^+ down to very 
small angular scales. By an unfortunate conspiracy (or maybe 
not) ^- falls off rather slowly towards large scales, and ^+ is 
dominant for ?? — > 0. Since a small-scale cutoff ^min due to 
close projected galaxy pairs, as well as a maximum scale ^max 
due to the finite observed field are unavoidable, an exact E-/B- 
mode separation using the above-defined correlation functions 
is not possible. 
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Fig. 4. The E-/B-mode con-elation functions. Top panels: 
Mixing between ^e+ and ^b+ due to a finite field. The curves 
correspond to an infinite field and to fields with maximal scales 
of 0max - 1 deg, 30, and 18 arc minutes as indicated in the 
panel. (The E- and B-mode for ^max - 6 deg is not shown 
in the upper panel.) Bottom panels: Mixing of ^e,b- due to a 
small-scale cutoff O^ia- In the small panels, the ratio between 
the corrupted and the true E-mode is shown. 



Similar to the case of the aperture-mass dispersion 
(Sect. |2li, we define the functions ^e,b+(0, ^max) and 
^E,B+(&, Ornm)- In the first case, the infinite integral for ^u 
JlOl l is truncated at ffmaxJ in the second case, the lower 
integration limit for ^l is replaced by O^ia- 

In Fig.|3] the mixing of the E- and B-mode correlation func- 
tions is displayed. The large-scale cutoff leads to an underesti- 
mation of the true E-mode on all scales, even down to arc sec- 
onds. At the same time, a spurious B-mode ^b+ is introduced. 
Even for a maximum angular separation of 1 .5 degrees, the E- 
mode is off by more than one percent on all relevant scales. The 
'minus' -E-mode ^e- suffers badly from the close pair cutoff on 
scales as large as several arc minutes. 

Figures|5]and|6lshow the deviation of the corrupted E-mode 
^E-^(^, ^max) and ^E-(S, ^min) from the true E-mode ^^+{ff) and 
^E-(0, respectively. 
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Fig. 5. On scales 6 smaller than 6^, the ratio ^e+(0, 0\ 

is lower than 1 - p, where lines corresponding to p 

and 1 % are plotted as a function of the maximum separation 

^max- Thick curves correspond to a mean source redshift ofz — 

1, thin lines toz = 1.5. 
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Fig. 6. On scales 6 larger than Oo the ratio ^e-(^, dmax)/^E-{0) is 
smaller than 1 - p, where lines corresponding to p - 10,5,2 
and 1 % are plotted as a function of 0niin- 

4. Dependence on cosmology and source galaxy 
redshift 

The amplitude and the shape of the shear statistics considered 
here depend, of course, on the underlying cosmology and the 
resdhift distribution of source galaxies. However, the ratio of 
the observed E-mode in the presence of the small-scale cutoff 
Omin to the true E-mode ((M^ ), ^e,b-) is very insentive to the 
convergence power spectrum. In the case of ^e.b+j the differ- 
ence between the E-mode corresponding to a maximal scale 
^max and the true E-mode depends only little on the underlying 
model. However, the amplitude of ^e+ depends quite strongly 
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on the amplitude of the power spectrum, which is governed in 
particular by Q.^, cr<i, and the source redshift distribution. As a 
consequence, the ratio of the observed to true E-mode depends 
on those parameters. For a deeper survey, this ratio will be 
lower than for a more shallow one. For all plots in this paper, if 
not indicated otherwise, a mean redshift of z = 1 .0 is assumed. 
The cosmology is ACDM with n^ - 0.3 and o-g = 0.85. 

5. Summary and conclusions 

We have quantified the mixing between the E- and the B-modes 
which arises due to the lack of information about the shear cor- 
relation on both very small and very large scales. The former 
limit 0min is due to close projected pairs, the shape of which 
cannot be determined reliably. The latter cutoff 0,nax is related 
to the finite size of the observed galaxy fields. Apart from that, a 
fundamental limit exists since only a 4;7r-survey can in principle 
yield an unambiguous E- and B-mode decomposition which, 
however, is unrealis tic due to the occlusion of parts of the sky 
by the Milky Way JBunri 120031) . We used the aperture-mass 
dispersion {M~ ^) and the shear-correlation functions ^e,b+ as 
measures of the E- and B-modes, respectively. 

Even if the cutoff due to close galaxy pairs is as small as 5", 
which is feasible for ground-based observations, the E-mode is 
underestimated by 10% (1%) on scales below V (3'4), using 
{M~) with the polynomial filter function. Moreover, a negative 
B-mode signal appears. Since even with space-based measure- 
ments, a cutoff of about 0,nin - 1" exists due to blending galaxy 
images, it is fundamentally not possible to know whether there 
is a B-mode present on scales smaller than one or two arc min- 
utes with very high (sub-percent) precision. This cutoff does 
not introduce a B-mode into the 2PCF, however, it prevents the 
clear detection of any B-mode still present in the data that may 
corrupt the 2PCF on much larger scales than 0^1^- 

The E-/B-mode correlation function ^e,b- behaves in a sim- 
ilar way than (M^ ). The '+'-version, ^e,b+i on the other hand is 
affected by a large-scale cutoff 0^^^ that causes a leakage from 
the E- to the B-mode on all angular scales. For fields smaller 
than 1 .5 degrees, the precision of the E-mode correlation func- 
tion is never better than 1%. 

In order to obtain knowledge about the B-mode in the data, 
one has to combine ^e,b+ and ^e,b-- For ground-based data 
(^min = 5"), ^E,B- allows One to constrain the B-mode at the 
percent level on scales larger than 3'3. To get this precision on 
larger scales using ^e,b+, the measurement of the shear corre- 
lation up to about 0„ax = 7 degrees is necessary. From space- 
based observations with flniin - 1"- ^max - 3.4 deg (0max - 2.4 
deg) is needed for a mean redshift of I = 1 .0 (z = 1.5), respec- 
tively. 

Up to now, various strategies have be en employed t o esti- 

I 1 J Tl 

mate the E- and B-mode from shear data. lMassev et al.lll2005ft 
calculate the E-/B-mode correlation functions by extrapolating 
^- beyond measured scales, using a theoretical model for Pk- 
However, this method explicitely makes the assumption that 
no B-mode is present on larger scales than probed by the sur- 
vey. Choosing a wrong model can change the amplitude of f e 
on the measured angular range, which might bias the cosmo- 
logical interpretation of the shear correlation. In particular, fini 



and o-% might be biased, which, although withi n the systematic 
errors for current measurements (see Fig. 10 o nHevmans et alJ 
l2005l) . will b e problematic for future sur veys. 

Recentlv. lvan Waerbeke et al.l ( 120051) use the aperture-mass 
dispersion to calibrate the E-mode shear-correlation function. 
Since their measured B-mode is consistent with zero on large 
scales (10' - 50'), they shift ^e,b such that the B-mode vanishes 
on corresponding angular scales (2' - 10'). It has to be noted 
that a B-mode on very large scales does not have any influence 
on (Mg ) on smaller scales, but it will corrupt the correlation 
function. Moreover, there is no clear correspondence between 
the angular scales probed by the correlation function and the 
aperture-mass dispersion. 

Alternatively, the aperture-mass statisti cs can be measure d 
and fitted with a cosmological model (e.g. ljarvis et alJl2003l) . 
In that case, either scales below 0q have to be omitted from 
the analysis, or the theoretical prediction of {M~ ) has to be 
obtained using eq. (|6jl with the cutoff fl^in specified by the ob- 
servations. Note, however, that in this case the B-mode signal 
should also be fitted with the predicted, non-zero {M±(0, 0mm)} 
if no lensing information is to be lost. 

At present, the prediction on the shear signal at angular sep- 
arations on arc-minute scales and below is very difficult and 
inaccurate. These small scales should not be disregarded since 
the shear signal on these scales is very high and contains unique 
information about the non-linear and non-Gaussian features of 
the large-scale structure. In the future, when our non-linear 
models of structure formation improve and baryonic effects are 
taken into account, shear measurements on small scales will be 
of great use. Since the shear signal on small scales is partic- 
ularly susceptible to contamination by intrinsic galaxy align- 
ment, a clear separation between the E- and B-mode on these 
scales is mandatory. Only then can one make certain that the 
measurement is free of systematic errors and has the necessary 
quality to be used for precision cosmology. 
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